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Black Hole Topology 



Investigations of classical signature change have generally envis- 
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Q\ \ aged applications to cosmological models, usually a Friedmann- 

Lemaitre-Robertson- Walker model. The purpose has been to avoid 
the inevitable singularity of models with purely Lorentzian signa- 
i . ture, replacing the neighbourhood of the big bang with an ini- 

i 

(3JT). tial, singularity free region of Euclidean signture, and a signature 

change. We here show that signature change can also avoid the 
; singularity of gravitational collapse. We investigate the process of 

j_j ■ re-birth of Schwarzschild type black holes, modelling it as a dou- 

ble signature change, joining two universes of Lorentzian signature 
through a Euclidean region which provides a 'bounce'. We show 
that this process is viable both with and without matter present, 
but realistic models — which have the signature change surfaces 
hidden inside the horizons — require non-zero density. In fact the 
most realistic models are those that start as a finite cloud of col- 
lapsing matter, surrounded by vacuum. We consider how geodesies 
may be matched across a signature change surface, and conclude 
that the particle 'masses' must jump in value. This scenario may 
be relevant to Smolin's recent proposal that a form of natural se- 
lection operates on the level of universes, which favours the type of 
universe we live in. 
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1. Introduction 



Space-time in general relativity is usually considered to possess a metric of Lorentzian 
signature. Positive definite metrics, with a Euclidean signature, have come into prominence lately 
through the Hartle and Hawking program concerning the wave function of the universe [1] - - [5]. A 
general aim of that program is to try get a handle on the boundary conditions of the universe, and 
an intriguing suggestion made in [6, 1] is that the universe has no boundary, i.e. no origin where 
initial conditions have to be set, which is only possible if space-time emerged from a Euclidean 
region preceeding a change of signature. Another interesting development is the introduction of 
Euclidean wormholes. These wormholes can arise in one universe and connect it either to itself 
or to another universe. In order to attribute a transition probability, for example, between two 
Lorentzian regions, integration of the action along the tube connecting the two regions under study 
is required. In normal Lorentzian space the path integral approach leads to oscillating behavior, 
and hence to non-convergence of the integral. To obtain convergence, the substitution t ^ it is 
applied, thus introducing a Euclidean signature. This means in effect that we have two Lorentzian 
regions connected through a classically forbidden Euclidean region. 

Paralleling the Quantum cosmology program, papers [7] - - [10] pointed out that the possi- 
bility of a change in the signature of the metric is not restricted to a quantum description of General 
Relativity. It was shown in [8, 9] that classical General Relativity does not prevent the existence of 
Euclidean regions and some examples of signature change in the Friedmann-Lemaitre- Robertson- 
Walker metric were produced. Further investigations were pursued in [11] - [30]. Even though 
the metric signature is invisible to the Einstein Field Equations, it should be noted [8, 19, 23] that 
a change of signature is not, either goo goes through zero, in which case the metric is degenerate 
and g 00 singular there, or both g 00 and g 00 jump from positive to negative values, in which case 
the metric is discontinuous. In either case, the Einstein Field equations cannot be defined in the 
usual way at the signature change [23, 27]. 

The Friedmann-Lemaitre-Robertson- Walker model has been considered repeatedly in both 
classical and quantum cosmological signature change, but few other models have been considered, 
especially in the classical signature change literature. This paper constructs a classical model 
of signature change within the black hole topology, using Schwarzschild and Lemaitre-Tolman 
models — i.e. a Kruskal-Szekeres type of manifold. It examines the transition from a black hole, 
through a signature change to a Euclidean region which reverses the collapse process, leading to 
a second signature change, and the birth of a white hole and a new universe. We also determines 
whether the signature change surface can be hidden inside the horizon. It continues the approach 
of papers [8, 9, 19] by exploring strictly classical signature changes in the Schwarzschild and 
Lemaitre-Tolman metrics. 

The investigation of transitions between Lorentzian geometries through a Euclidean region 
are also of interest when considered in conjunction with Smolin's idea [31]. Smolin's hypothesis 
is a proposed mechanism for determining the particular values of fundamental physical constants 
observed today, and thus justifying the anthropic principle. In Smolin's paper, life supporting 
characteristics are linked to the existence of stars whose abundance is linked to the abundance 
of black holes. It is envisaged that each 'universe' either expands and re-collapses or expands 
indefinitely, possibly forming one or more black holes. Instead of classical singularities occuring 
- a crunch singularity or black hole future singularities — quantum cosmological tunneling gives 
birth to new universes, and hence a 'natural history' of universes arises. Furthermore, Smolin 
proposes that the process of tunnelling generates small random changes — 'mutations' — in the 
values of the physical constants. Those combinations of values for which the universe generates 
many black holes will lead to large numbers of offspring having very similar values. Thus, after the 
passage of many generations of universes, the population of universes will come to be dominated by 
those that generate lots of stars and black holes. This parallels natural selection in that the 'fittest' 
universes reproduce prolifically, but differs in that all blood lines (sets of constants) survive. 



In what follows, we consider the junction conditions at a signature change in a Schwarzschild 
type metric, the choice of signature change surface, the form of the metric in the Euclidean region, 
and how geodesies should be propagated through. Generalisations lead naturally to the Lemaitre- 
Tolman metric and its Kantowski-Sachs limit, which allow more interesting results. 

The first step one has to take is to ensure that the various regions composing the space 
match geometrically. We adopt the Darmois junction conditions [32], and the application of them 
to signature change as presented in [19]. 

2. Junction Conditions and Conservation Laws 

Conventions 

We here work with 4-dimensional manifolds of Lorentzian and Euclidean signatures ( — h 
++) and (+ + ++) respectively. Greek indices range 0-3 and Latin indices 1-3. Subscripts E 
and L denote quantities defined in or evaluated in Euclidean and Lorentzian regions respectively, 
and expressions without such subscripts are valid in either region. These may also be written as 
superscripts, to avoid confusion with tensor indices. Geometric units are used, G = c = 1, and 
the cosmological constant is set to zero, A = 0. 

Darmois Matching conditions 

In standard Darmois matching, where no signature change occurs, space is composed of 
two regions, V + and V~ , with a common boundary surface E. More precisely, an isomorphism 
■if) : E + — > E~ allows us to identify the boundaries of V ± , E + = E~ = E . The two regions have 
coordinate charts x+ and xt and metrics g+ u and g~ respectively. Setting the intrinsic coordinates 
of the junction surface to be C+ = €- = the locus of the surface is given parametrically in V ± 
by x± = x±(£ l ), or by E ± (x±) = 0. We write Q ^ to denote evaluation of some quantity Q in the 
limit as the surface is approached from either region, and [Q] to denote the difference between the 
two limiting values 

[Q]=Q\+ -Q\- (1) 

The Darmois conditions [32] require the continuity of the first and second fundamental forms of 
the junction surface — i.e. the intrinsic metric and the extrinsic curvature. The intrinsic metric 
is obtained by projecting the 4-metric onto E using the basis vectors ef of E 



dx» 



lij — 9/J.v e i e j j e i — ati (2) 



The extrinsic curvature describes the surface's shape in the enveloping space, and is the projection 
onto E, of the rate of change of the surface normal n M in the enveloping space, with respect to 
position on E. 

K ll = (V l ^)e 1 e^-n^— + T^ ww j, n, = (3) 

where e n = n v n u = +1 if E is time-like and —1 if it is space-like. To conform with (1), the sign in 
(3) is set so that the n± point from V~ to V + on both sides of E. The Darmois conditions may 
now be given as: 

[h 3 ]=0 k [tfy]=0 (4) 

In the constant signature case, according to [33], these are equivalent to the Lichnerowicz matching 
conditions [34], whereas the O'Brien and Synge conditions [35] are too restrictive. 

When we introduce a signature change at E, the equivalence between the Darmois and 
Lichnerowicz conditions breaks down. Both the Lichnerowicz and O'Brien and Synge conditions 
insist that all the 4-d metric components be matched on either side of the junction surface, 
leading to a degenerate metric, a non-affine time coordinate, and breakdown of the Einstein field 



equations. We select the Darmois matching conditions as they are invariant to the coordinates 
chosen on either side. They require no modification at a surface of signature change. In fact 
they are blind to the change of signature, thus extending the signature blindness of the Einstein 
field equations. A signature change surface is necessarily space-like, so e n = n^n^ = +1 in the 
Euclidean region, and —1 in the Lorentzian region. 

Conservation Laws 

In [19, 20] the implications of signature change for conservation laws are worked out. Con- 
servation laws are based on the divergence theorem — i.e. the components version of Stokes 
theorem for a 3-form in a metric space. The theorem requires a region W, bounded by a closed 
surface S with outward pointing unit normal m a , smooth non-zero volume elements d A W and d 3 S 
on W and S respectively, a smooth non-zero metric, so that the inverse metric is well defined, and 
a smooth field ^/ s : 

I ¥m s d 3 S = [ V s V s d 4 W (5) 
Js Jw 

It should be noted that n a is the normal to the junction surface E, and m a is the normal to S, 
the closed boundary of W; 

These conditions are not satisfied through a signature change. Thus physical conservation 
laws need to be revised. For the electro- magnetic field we work with the 4-current = J s , and for 
the gravitational field, a component of the Einstein tensor ty s = G 7<s t> 7 where v 1 is some suitable 
smooth vector field. Since v 1 and v 7 are not both smooth through a signature change, = G^v 7 
is also considered. 

Firstly, at a boundary where no signature change occurs, the Darmois junction conditions 
may be used to patch together two regions that adjoin the boundary on either side, and within 
which the divergence theorem does hold. It is shown that these conditions, which give rise to 
Israel's identities [36] for the Einstein tensor, 

[G^rf] = (6) 

[GVef^l = (7) 

where 

G^n" = ^{K 2 - - e n &} (8) 

G^n v = V^f-V^ (9) 

¥? and being the intrinsic curvature invariant and covariant derivative of the 3-surface, and 
K = g^Kij = Kj, are sufficient to ensure conservation of energy-momentum through E. Combined 
with suitable junction conditions on the electro-magnetic field, they also ensure conservation of 
4-current, with similar results applying to other fields. 

At a surface of signature change, e n now flips sign across E, and this leads to modified 
Israel identities 

\G^n v \ = -ft (10) 

[G^rf] = (11) 

It is necessary to distinguish two normals to E: 1$ = /dx 5 = e° s and n s = dx s /<9£° = where 
n 7 n 7 — e n — l~p and g-^n 5 = e r X 1 . A similar analysis of the divergence theorem through E is 
made, paying careful attention to index position, the definition of the extrinsic curvature, and 
the directions of the various normals. It is found that, in the process of patching together the 
two divergence theorems on either side of the signature change, the combined theorem aquires a 
surface term, so conservation laws must in general be modified. The result is 

VPmpcPS-J Ed 3 S Q = J V^d 4 W (12) 



where 

E = (*« J+ - = [*X] (13) 

and S is the region of £ enclosed by S. For each of the four choices ty s = G 1& v 1) t> 7 = Z 7 or 
and \P 5 = G 7 f 7 , t> 7 = n 7 or e] the surface term E = E(vj) or £(t> 7 ) in the conservation law is 

-R (14) 

2(V j K ij - ^ ij3 VjK) (15) 

^ 2 - KijK ij (16) 

2(V^ - V,K) (17) 



E(y = [c^u^] = 

E(eL) = [G°%4] = 

£(n Q ) = [G^ Q r/] = 

E(e?) = [CXef ] = 



The main results are expressed in a way that allows any set of junction conditions to be 
applied at the signature change, so this conclusion is independent of choice of junction conditions, 
as well as being coordinate invariant. Alternative approaches to signature change, which emphasise 
maximum smoothness of the metric and the matter, are able to eliminate some, but not all of the 
surface effects. Removal of all surface effects requires that the surface of signature change not only 
have zero extrinsic curvature, but also have zero (3-d) Ricci scalar, which eliminates all realistic 
cosmological models. 

At a signature change then, the Darmois conditions still impose the same number of metric 
conditions as was sufficient for no signature change, and they still result in a modified set of 
conservation laws, albeit with surface effects. These can be understood as a consequence of the 
change in the character of physical laws. Whilst it is of interest to follow the maximum continuity, 
Lichnerowicz type approach, it is argued that physically interesting scenarios may be eliminated 
by it. The following models are an example. 

3. The Schwarzschild Case 

In the usual Schwarzschild line element [37] the signs of the metric components gxT and gRR 
interchange across R = 2M, thus leading to reinterpertaion of the roles of R and T. Consequently 
it is not clear which sign we should change to introduce the signature change, assuming that the 
general form of the metric is retained. We shall investigate both possibilities, so we insert two 
new sign factors €t — ±1 and e R = ±1 in the metric: 

ds 2 = e T (l- ^) dT 2 + e R (™- 1) _1 dR 2 + R 2 dQ 2 (18) 

where dVt 2 = <i# 2 +sin 2 9 dcf 2 . We are interested only in transitions from the standard Schwarzschild 
metric to a Euclidean region, so we disregard the sign combination which gives us a second (non- 
vacuum) Lorentzian manifold (e^ = +l,e# = +1). With both the other sign combinations 
(st = — 1, £r = +1 and 6t — +1, £r — — 1) it appears that there is a Euclidean region on one side 
of R = 2M, and a "double Lorentzian" (or "Kleinian") region (with two time-like coordinates) 
on the other side. It will become clear that the two Euclidean 'regions' are in fact geodesically 
complete manifolds. 

Whilst Euclidean regions have no time, there will be a direction which is the extension of 
the time direction in the Lorentzian region. One may determine whether the Euclidean metric is 
'static' or 'dynamic' relative to this direction. 

The Einstein tensor components for this metric are [38] 

_ -e T (l + e R )(l-2M/R) _ (1 + e R ) 

Gtt- ^ , G RR - R2(1 _ 2M/R y Gee-0, G„-0 (19) 



A vacuum solution requires = —1. Transitions requiring a change of sign in €r introduce non- 
vacuum solutions with strangely behaved matter, which we shall consider, since we don't really 



know what to expect in Euclidean signature physics. The Ricci scalar and the Kretschmann scalar 
are 

R£= 2 -^, k = R^R^ = ^Smi^MM±^l , 20 ) 

where R^ and R^\a are the Ricci and Riemann tensors. Regardless of signature, k is always 
singular at it! = 0, and never at R = 2M, so the Euclidean region (er = +1, e_R = —1, R > 2M) 
has no singularities. 

We now investigate whether signature change is possible on two simple spacelike surfaces 
in R > 2M a constant T surface, and in R < 2M a constant R surface. More general transition 
surfaces will be considered using a different metric form. The surface coordinates may then 
be chosen to be identically 3 of the enveloping coordinates in V~ and V + - - viz (T, 6, <fi) or 
(R : 9,(f)). (c.f. [39] in which the Euclidean solution for vacuum with a cosmological constant is 
found.) 

Constant T surface 

The intrinsic metric, unit normal, and extrinsic curvature of S are given in either region 

by 

do* = {1 _~2M/ R ) dR2 + R2dn2 > z T e n = l (21) 

am 

ra" = = = , n »n fM = e n = e T , K tj = (22) 

U T e n (l - 2M/R) 



and the surface effects are all zero. 

E(L) = 2{l R2 R) = °> £(4) = 0, % a ) = 0, £(e?)=0 (23) 

The two choices of future direction for n M are equivalent in a static metric. The choice of a 
standard Schwarzschild solution in the Lorentzian region sets the sign of €r to —1, so 6t flips 
across X, and this requires that Rl > 2Ml for a space-like surface. Although (21) is singular at 
R± = 2M±, all constant T surfaces on the Lorentzian side intersect that point, which is only a 
coordinate singularity, being the middle of the Schwarzschild wormhole at its moment of maximum 
expansion. Applying = and requiring the angular coordinates on either side to coincide, 
@e = @l, <Pe = also fixes the areal radius and mass terms to be the same, Re = Rl and 
Me = Ml- Obviously [Kij] = imposes no further constraints. Since all T = constant surfaces 
are equivalent (for a static metric), this result is not surprising. 

This matching corresponds to vacuum both before and after the signature change, but the 
change surface extends to spatial infinity in both exterior regions Rl > 2Ml- We are really seeking 
a change surface that is near the singularity R = and hidden inside the horizon. No spacelike 
surface is further from R = or less hidden than that of constant T. 

Since the middle of the throat at maximum expansion is a stationary point of the Killing 
vector x 11 = $t on the Lorentzian side, and since all Lorentzian constant T surfaces match to all 
Euclidean constant T surfaces, this must be a stationary point on the Euclidean side also. This 
leads us to suspect that it is not possible to find two separate, non-intersecting T = constant 
surfaces in the Euclidean region. In other words, we can't construct a Euclidean region between 
two separate Lorentzian regions. 

Constant R surface 

This is the simplest non-vacuum case. The fundamental forms and surface effects are: 



da 2 = e T (l - 2M/R)dT 2 + R 2 dQ 2 , e R e n = 1 



(24) 



n" = -^/e R e n (2M/R-l)5^ n"n, = e n = e R (25) 



-e T MJe R e n (2M/R-l) 



R 2 ' 

K ee = -R^Je R e n (2M/R-l), = sin 2 9 K 9d (26) 

E{Q = ^, E&) = 0, E(n a ) = ^^ = ^, E(et) = (27) 

where we have chosen n M to point in the direction of collapse, i.e. towards R decreasing. A similar 
analysis gives us = —1, Rl < 2M L , and e R flips; choosing 9 E = 9 L , (j) E = 4>l and T E = T L with 
fyij] = =>- Re = Rl, and M E = M L . No further restrictions are necessary to ensure [K^] = 0. 

This demonstrates that matching can be achieved on a surface that is entirely inside the 
horizon, but since the Euclidean region is not empty in this case, it still leaves open the in- 
terpretation of the energy stress tensor on the Euclidean side. Further, since R is the timelike 
coordinate on the Lorentzian side, R is the nominal 'time' direction on the Euclidean side too, 
being orthogonal to the transition surface, so this Euclidean manifold, with R < 2M, is 'dynamic'. 

Geodesic Coverage 

What do the particle paths look like in the combined space? We now investigate the 
behavior of radial timelike geodesies. Adding the angular components of the motion should present 
no problem, as 9± and 0± are identified at X. One aim is to verify that the space is geodesically 
complete, and in the context of this paper geodesies that end on a curvature singularity are 
considered to be as complete as is possible. The second aim is to see how geodesies should be 
continued at the transition, and whether the set of all geodesies arriving at the Lorentzian side 
of the transition generate all possible geodesies emerging on the Euclidean side. Three schemes 
for continuing geodesies are considered. Two of them attempt to match particle 4-velocities (unit 
normal tangent vectors), and one attempts to match 4- momenta. 

The geodesic equation 

M ^VX = (28) 
with the condition for a 'time-like' unit normal 

= e « ( 29 ) 
where e n = at a constant T transition or e R at a constant R transition, leads to the acceleration 

and gives the following unit tangent vectors, where the signs have all been chosen so that positive 
h and q values always give consistently future directed infalling tangent vectors where T > 0. 

h, 



Lorentzian: ^ = { {1 _ 2M/R) , ~^h\ - (1 - 2M/R) , 0, Oj , q L = ±1 (31) 

There are three types of geodesies: (a) 1 — 2M/R < h\ < 1: geodesies recollapsing from past 
to future singularities, R = 0, with a maximum at R = 2M/(1 — h\) > 2M; (b) h\ > 1: 
monotonically ingoing or outgoing geodesies with finite velocity \Jh? L — 1 at R — oo, reaching 
R = either in the past or future; (c) h\ = 1: marginal monotonic geodesies with zero velocity 
at R = oo. 

h f 



Euclidean, R > 2M: u» = ^ 1 _ 2M/R) , ~Qe y/(l - 2M/R) - h% , 0, Oj , q E = ±1 (32) 



In this case there is only one type of geodesic, descending from R = oo, through a minimum 
at R = 2M/(1 — h 2 E ) > 2M, and re-expanding back out. The allowed range of He is then 
< h 2 E < (1 — 2M/R), and all geodesic paths are restricted to the region R > 2M. The only 
geodesic reaching R = 2M is the one with He = which in effect is a stationary point. This is in 
accord with the range of R at E, and confirms that the region R > 2M is a geodesically complete 
manifold. 

Euclidean, R < 2M: u» = ( 2M ^_ 1 , -q E \/(2M/R - 1) - h\ , 0, 0^ , q E = ±1 (33) 

These geodesies all expand from R — and re-collapse back, and have maxima at R — 2M/(1 + 
h E ) < 2M, where h E < (2M/R — 1). Thus the region R < 2M is also geodesically complete, but 
does not have the desired bouncing property. (R = 2M would not actually be encountered in this 
region of our model, since R% < 2M.) 

Since geodesic tangent vectors have a unit magnitude that flips sign across a signature 
change, it is impossible to match all components of u± across E. We consider here three possible 
schemes, and summarise the resulting conditions in Table 1 below: 

(i) Match u R values, with u^u^ = e n giving a jump in the u T values; 

(ii) Match u T values, with u^u^ = e n giving a jump in the u R values; 

(iii) Match all components of the 4- momentum P M = mu^, allowing P^ and 
\P^P^\ = m 2 to jump. (This is equivalent to matching non-normalised tan- 
gent vectors, for which the metric degeneracy at E is irrelevant, and standard 
existence and uniqueness theorems guarantee that geodesic continuation of all 
tangent vectors is possible.) 



TABLE 1 GOES HERE 



TABLE 2 GOES HERE 



When these are compared with the allowed ranges of h in each region, summarised in Table 
2, we see that conditions (i) and (ii) do not allow the continuation of all possible geodesies that 
might arrive at either type of signature change surface, whereas (iii) continues all geodesies at 
both types. For example, identifying u E — u R at a constant T surface gives us 

h 2 E = 2(1 - 2M/R) - h 2 L , R>2M, h\ > (1 - 2M/R), h 2 E <(l-2M/R) (34) 

so we can easily find a large enough value for h\ to make h E imaginary. For case (iii) we find 

Constant T transition: ( ^ V = ( M = ^ - - - 1 (35) 

\m L J \h E ) (1-2M/R) V ; 



Constant R transition: ( — — ( t^-\ = , , ttM r + 1 (36) 

\m L ) \h E j (2M/R-1) V ; 

Since h\> {1 — 2M/R) the particle's Euclidean mass is always greater than it's Lorentzian mass 
rriE > rriL, as well as hi > h 2 E . 

Thus it turns out that condition (29), forcing the tangent vectors to be unit vectors, is too 
strong an assumption, and doesn't permit all particle paths to be matched through E. Rather, 
the matching of geodesic 4-momenta is the only way of extending all particle paths through a 



signature change. The conclusion that the 'rest-mass' parameter of a particle has to jump, is 
consistent with the fact [19] that in general tensors cannot be smooth through E in both covariant 
and contravariant forms, and that the density can jump across a signature change. In Fig. 1 we 
summarise the properties of geodesies arriving at £ at some particular value of The horizontal 
axis is the parameter h\ and the plot covers a representative range of permissible h\ values for 
the chosen 



FIG 1 GOES HERE 



Summary 

Within the Schwarzschild metric form, signature change is possible on constant R surfaces 
inside the horizon, but the resulting Euclidean region has strange matter, and continues to collapse 
to a singularity. A second signature change back to a Lorentzian region is of course possible, but 
only at a smaller R, closer to the future singularity. Signature change is also possible on constant 
T surfaces, leading into a Euclidean region which is vacuum and has geodesies which do bounce. 
But constant T surfaces are entirely outside the horizon, and so not very interesting, since a second 
transition back to a Lorenzian spacetime results in the same future as the one that was avoided. 

Attempts to continue geodesies through a signature change indicated that one must match 
the 4- momenta, which means that particle rest masses have to jump. If a second signature change 
back to a Lorentzian metric occurred, the particle mass would not return to its original mass 
unless the model were highly symmetric. 

4. The Lemaitre-Tolman Case 

We now shift our attention to the Lemaitre-Tolman metric [40, 41], primarily because it 
allows us to deal simply with more general surfaces in spherical vacuum, and secondly because it 
makes possible a generalisation of the black hole topology to non-empty models, thus describing 
more realistically the collapse of matter into black holes, as well as the more standard cosmological 
collapse of matter, where no wormhole topology is involved [42]. This gets us closer to Smolin's 
scenario. 

In the vacuum case, the Lemaitre-Tolman metric with appropriate choice of parameters 
can describe the full Schwarzschild-Kruskal-Szekeres manifold, avoiding a coordinate singularity 
at R = 2M and the accompanying change of character of the Schwarzschild R and T coordinates, 
and making it clear which metric element should change sign at a change of signature. Its two 
arbitrary functions make it much more flexible than the Kruskal-Szekeres metric. 

The diagonal, synchronous, spherically symmetric metric, with an added factor of e = ±1, 

ds 2 = edt 2 + B 2 (t, r)dr 2 + R 2 (t, r)dVl 2 (37) 
leads to the following Einstein tensor: 

e(2BRR" + BR' 2 - 2B'RR' - B 3 ) + (2B 2 BRR + B 3 R 2 ) 

Gtt ~ w& (38) 

G, = (39) 
(R' 2 -B 2 ) + e(2B 2 RR + B 2 R 2 ) 

Gtt - ^ (4UJ 

G H _ R[(BR" — B'R') + t(B 2 BR + B 2 BR + B 3 R)] 

^e = Gee = P (41) 

(42) 



where ' = d/dr & '= d/dt and the cosmological constant is taken to be zero. 

Solving the Einstein field equations for co-moving matter, w M = <5f , and zero pressure, 
p — 0, gives the Lemaitre-Tolman model, and we get 

B 2 = /(r)>-l (43) 

-ete = 2 -f + f{r) (44) 
M 

R = e- 2 (45) 

2M' , x 

87rp = G tt = -e— (46) 



, . 2B'i2' (i?') 2 2i?" \ (1 - eR 2 ) /2i2 2BR B 
2M' 



_B 3 i? ^i? 2 E 2 /? / i? 2 \R BR BJJ 

(47) 




2\ 



2 (J^ BW_ eBRY (R^_ BR' \ 
+ \B 2 R B 3 R + BR ) + € \BR B 2 RJ 

/3M /2 8M'M 12M 2 \ 

\R^~l¥R T + ~R ir ) ( ' 

where / = f(r) and M = M(r) are arbitrary functions of coordinate radius r, p is the density, 
and k is the Kretschmann scalar. Singularitites in k and p occur at it! = and R' — regardless 
of e. Shell crossings occur where i?' = 0, since shells of matter at a different constant r, arrive at 
the same areal radius R(t, r) and intersect each other. In vacuum, M' = 0, so p is zero and k is 
finite, and there is no physical problem, but there is a bad coordinate coverage of the space. In 
non-vacuum cases care needs to be taken to select the arbitrary functions which do not give rise 
to these physically troublesome shell crossings [43]. 

In the standard Lorentzian case (e = —1), f(r) is a kind of local energy constant which 
determines the type of time evolution — elliptic, parabolic or hyperbolic — as well as the local 
geometry, and M (r) represents the total effective gravitational mass within comoving radius r. 

In the Euclidean case (e = +1) the acceleration (45) is everywhere positive, provided we 
select a positive 'mass' term, so a bouncing Euclidean universe can be achieved. This requires 
to be negative in order to keep R real, i.e. — 1 < < 0. 

We now obtain solutions to the evolution equation (44), in terms of a parameter r], and 
a = a(r), a third arbitrary function of r, which is the time of the big bang R = 0, or if we use the 
time reverse of the following equations, the time of the big crunch. Solutions with / > or M < 
are discarded. Although we can't be sure negative 'mass' solutions are physically disallowed in a 
Euclidean manifold, they all reach the R = singularity, and none of them bounce (re-expand), 
so they do not serve our purpose. 

Lorentzian region e = — 1: elliptic solution, — 1 < f L < 

M L M l 
R ( t ' r ) = 7~~ TTV 1 -cosr] L ), t= (r] L - sin r] L ) + a L (r) (49) 



Euclidean region e = +1: M E > 0, -1 < f E < 

M E M e 
R(t,r)= ^ (coshrjE + 1), t = ^ 3/2 (smh r] E + r] E ) + a E (r) (50) 

Any Lemaitre-Tolman model with M' = is a vacuum model, and thus for e = — 1 repre- 
sents at least a section of the Kruskal-Szekeres-Schwarzschild space time in geodesic coordinates. 
However not every selection of the arbitrary functions gives complete coverage of the manifold. 
Novikov coordinates [44] do cover the entire manifold, and are obtained with the following choices 

Ml = const, f = ^ - ( '^ )2 , a L (r) = (51) 

for which the surface t = is a simultaneous time of maximum expansion, and /(0) = —1 at 
the Schwarzschild throat, increasing monotonically to as r — > ±oo. This topology — two 
sheets joined by a throat — may easily be extended to non-vacuum everywhere [42] by setting 
Ml = Ml(t) with a minimum value at the throat. It is the form of f(r) which determines the 
topology. If the asymptotic regions are closed FLRW cosmologies (/ = —kr 2 , k = +1), then we 
still expect / to rise very close to zero before decreasing again. In such dense black holes, the past 
and future event horizons are split, and R = 2M L is an apparent horizon [42]. 

Matching conditions 

We perform the matching on the simplest possible surface, that of constant time, t = 
constant. In vacuum this is merely a coordinate restriction and not a physical one, because the 
origin of the time coordinate, a (r), is an arbitrary function of position. It amounts to finding the 
family of geodesies orthogonal to the transition surface, and using these as lines of constant r. 
The intrinsic metric of such a surface is correspondingly simple: 

da 2 = B 2 dr 2 + R 2 dtf = dr 2 + R 2 dQ 2 (52) 

1 + / 

When matching, a reasonable choice is to equate the angular parts, and to re-scale the coordinate 
radii, so that 

E = L , (p E = 0l, r E = r L (53) 

and = fixes 

Re = Rl = Rz, Be = Bl = By, (54) 

Since R is continuous across the junction and is a function of r only on S, i.e. R% = Rj;(r), we 
have also that 

R' E = R' L = R'z, =>f E = f L = f (55) 
Because of (43) the normal and the non-zero elements of the extrinsic curvature are: 

= y/ee^5p, n^n^ = e n = e, (56) 

K rr — y/ee^BB — 7-7—7, K ee = ^T n RR = RR = ^- (57) 

1 + / sin 



and [Kij] = leads to 
The surface effects are 



Rl = Re, Bl = Be R'l — R'e (58) 



, 7 ^ ^f(R') 2 2B'R> 2R" 1\ // /' \ 

B 2 E(K) = E{e«) = 4v /^ - ^) (60) 

E(n°) = ^ (61) 



The principal feature we are looking for is a bouncing universe, meaning a Lorentzian region 
matched to a bouncing Euclidean region that in turn may be matched to another Lorentzian 
region. This involves establishing the existence of at least two solution surfaces in the Euclidean 
region of the model under investigation. In general, given two space-like hypersurfaces, there will 
not be any geodesies that are orthogonal to both, so requiring both to be t = constant surfaces 
in the same coordinate system could well be restrictive. 

General transitions 

Five arbitrary functions, f(r), M E (r), M L (r), ax(r) and a E (r), are as yet unspecified. We 
now derive the necessary relations between them at a surface of signature change. We do not 
assume vacuum at this stage. Only models with / < and M E > give rise to a Euclidean region 
with a bounce. 

Condition (58) for R 2 combines with the evolution equation (44) to give 

M L -M E Ililh: or M L + M E = -fR E (62) 
The sign of R must still be matched. 

Inserting the parametric expressions for i? s (49) and (50) into (62) gives 

cos riL^ = ~j^, coshes = — (63) 

M E <M L , -1 < cos?] L <0 , 1 < cosh r] E < oo (64) 

which, combined with the continuity of R for collapsing models, yields 

*e < : %s = 2vr - oosT\-M E /M L ), r} E11 = - cosh" 1 (M L /M E ) (65) 

Since the Euclidean region doesn't have arbitrarily large R, the transition cannot happen arbi- 
trarily close to R = 0. On the transition surface t% is constant, so t)ly, and rj E z are functions of r 
only. Thus 

M E 

tET, = t-E = a E (r) + ^_^ 3/2 ( smn Vet, + Vet) (66) 

M L 

£le = *e = a L (r) + (Vls - sin?? LS ) (67) 

If possible, we choose a E (r) = 0, so that (66) gives a time-symmetric coordinate coverage in the 
Euclidean region. This permits a second copy of any transition surface found away from t — 0, 
and thus ensures a bounce. To obtain a specific solution, we fix ts and any two of Ml, M e , a E , a E , 
to obtain the others. 

Further, (62) plus the requirement that the transition surface be inside the external horizon 
R = 2ML m ax, where M^ max is the total exterior mass of the collapsing cloud, gives us 

_M L (r) + M E{ r) _ 

Zjlvl Lmax Zjlvl Lmax 

since < M E < Ml- At the centre we have /(0) < — {Ml m in/ ^ Ml max) and outside the cloud or 
at large r, f(r > r max ) < —(1/2), where r max is the smallest radius for which M^(r) = MLmax- 
This is a very stong restriction on /. 

Vacuum to Vacuum 

We set the mass term constant (hence p — 0) and equal on either side M E = M E . By (62) 
this gives i? s = 0, and by (49) and (50) the areal radius can only be matched at 

^e = 0, 77 LS = 7r, J R S = ^-, -l</<0 (69) 



the loci of minimum and maximum expansion of the Euclidean and Lorentzian coordinates re- 
spectively. Clearly, this case is equivalent to a constant T transition in the Schwarzschild case, as 
£ touches R = 2M but otherwise lies entirely outside the horizon. If we set a E = 0, to obtain 
a symmetric coverage of the Euclidean region, we have that the transition time is t% = — i.e. 
there is no 'time' between the two transitions. This confirms our earlier suspicion that mini- 
mum expansion at the middle of the throat is also a unique event in the Euclidean Schwarzschild 
topology. 

Retaining M E and M L constant, but not necessarily equal, (63) shows that both rj E and tjl 
are constant on the transition surface. The parametric expressions for t (66)-(67) then establish 
the relation between f(r) and a(r). Again symmetric coverage of the Euclidean region — a E = 
in (66) — would require / = constant and hence R% = constant. The only way to get surfaces 
which have constant t, R, and /, is to set / = — 1 — dealt with next. 

Constant R 

The constant t surfaces can also be made constant R surfaces, thus yielding the closed 
Kantowski-Sachs model [45] in Lemaitre-Tolman coordinates [46]. This is done by setting 

M = Mi j y/l + f dr + M , a = a ± J y/l + f dr + a , M , M u a , a ± constants (70) 
and then taking the limit / — > — 1, leading to 

Sill 71 t 

e = -l: R L = Mo L (l-cosr] L ), B L = 2M 1L - (M 1lVl + a lL )- ^— (72) 

(1 - COST] L ) 
Slllll 77 fp 

e = +l: R E = M os (cosh?7 £ + 1), B E = 2M 1E - (M 1E r) E + a 1E ) — '—— (73) 

(cosh?7 E + 1) 

In Lorentzian vacuum, M 1 = 0, e = —1, these coordinates only cover R < 2M^, and may be 
similarly incomplete in dense models. Shell crossings may be avoided for — 2n < (an,/ Mil) < in 
Lorentzian regions, but not in Euclidean regions. However, for a\E = 0, shell crossings are removed 
if transitions happen at | t] E t, \< 2.3994 which is the positive root of 2(cosh 77^ + 1)— r] E mih.r] E = 0. 

By (62), Rx = M QE + M 0L and for M 0E > and £ hidden (i? s < 2M 0L ) we need M 0L < 
Rz < 2M 0L < M 0E < M 0L . Equations (62)-(63) and (66)-(67) become 

Rx = M E + M 0L , cost ?ls = -^, cosh^ = ^ (74) 
a E + M 0E (smh.r]ET, + Vet,) = *s = a 0L + M 0L (i] LTl - shirks) (75) 

and because the matching of Ry, and R% no longer ensures B% and B% is matched, we get two 
extra conditions 

(M 0E M 1L - M 0L M 1e )Rt+ 
((Mil - Mib)* s - a 0L M 1L + a 0E M 1E + a 1L M 0L - a 1E M 0E )^M^ L - M$ E = (76) 
(Mil + M 1E )t E - a 0L M 1L - a 0E M 1E + a lL M 0L + a 1E M 0E ) = (77) 

Simplifications are obtained by requiring a symmetric Euclidean region a^ E = 0, a± E = 0. Vacuum 
to vacuum is not possible since M 1L = 0, M 1E = implies a\ E = = Oil, and thus B L = = B E 
at all times. Similarly, vacuum to non-vacuum is not possible. The dense models are highly 
symmetric, as they have uniform density on constant R surfaces. A sample set of values are: 
M 0E /M 0L = 0.1, 77 £S = -2.9932, r/ L s = 4.6122, a 0E = 0, t s /M 0L = -1.2943, a 0L /M 0L = 
-6.9015, Rx/M 0L = 1.1, a 1E = 0, M 1E /M 1L = 0.0407, a 1L /M 1L = -5.5476, 5 S /M 1L = 
1.1539, p L Ml L = 0.0570, p E Mh = 0.0057. 



Dust to Dust — Wormhole Topology 



We set Ml = M^(r) and Me = M^(r). Since space-time is no longer empty, not all geodesic 
coordinate systems are equivalent to the comoving one, so a symmetric coordinate coverage of the 
Euclidean region becomes essential to ensure the existence of a second transition surface. From 
(63) and (66) with a E = we obtain: 

jp ( \ _ (-/) 3/2 ^s smiles + t]ej: _ . , . . 

F{r) = — = = D(r]Es) 78 

M L cosh^ s 

where the right hand side defines D(t] E t,), and the left hand side defines F(r). For a Kruskal- 
Szekeres-Schwarzschild type topology (Lorentzian) [42], we expect / to take a Novikov-like form, 
i.e. symmetric, /(— r) = f(r), with /(0) = — 1 a minimum at r = 0, and rising monotonically. To 
cover the asymptotically flat regions at large R requires /(±oo) = (since f E > doesn't give a 
bounce). We expect the mass Ml to be minimum at r = and rising monotonically to a finite 
value. For example 

J = -7-, — 9 ' Ml = r - ; — 2 v y J 

1 + r 2 1 + r l 

With these choices and M L min = (l/3)M ima;E , F(r) and D(t)et,) are plotted in Figs 2a and 2b. 
The main features of F(r) are dictated by the topology, and are independent of the particular 
choices of f(r) and Ml(t). 



FIGS 2a AND 2b GO HERE 



The mapping between F and D is needed to fix the r dependence of r]Es(r). We have at our 
disposal only one constant, t^/MLmin, which is freely adjustable, so we need to select the section 
of the D graph which includes 0, in order to accomodate / — > 0. Also since F is monotonic, the 
range of tjes cannot extend through the maximum of D(rjET,), D max . Hence we have a restriction 
on when in the Euclidean evolution the transition can occur. To obtain /(0) = -1 we need 

/ = -1 =>- ts = sinh VEX + VEi; , go v 
M L coshes 

and consequently the following ranges are allowed: 

0<Y^<D max , < 77 £S < r] E1 : max , — <M E <M L (81) 

M L min COShrjEHmax 

where 

D m ax = 1-5434, at r] ESmax = 1.5434 = D max (82) 

Although the range f]ET,max < Vet, < 00 is not obviously precluded in principle, it results in a 
different sign for dr]Ej:/dr, which affects R^, pet, and plt, as well as limiting the range of R. 

We now find some sample solutions numerically. Our plotting procedure is as follows: 

— Select the Lorentzian mass, M L ] 

— Choose a transition time t% which complies with (81); 

- Generate values of M E which span all values allowed by (81), for the given t s and M L ; 

— For each Me value calculate: 

- / from (66): / = -[(sinh^ + ^ E )M s /(t s - a E )] 2/s ; 
r from (79): r = yj(l + /)/(-/); 

- i? E from (62): i? s = (M L + M E )/(-f); 

a L from (67): a L = t s - (r) L z - sinr7 LS )M L /(-/) 3 / 2 ; 



- Pex and p LS from (46): p s = 2M' jR^R^ 
(The last requires the values of M' E and R' , obtained from the derivatives with respect to r of 
(66), (67), (63), and (62)). 

Having found two surfaces where a signature change could occur, i.e. t^ E = — tfp L > 0, 
the idea is to excise the future singularity in one Lorentzian region, and the past singularity in 
the other, and join them with the Euclidean region. 

The following three models are typical. They use the forms Ml = {ML m i n + ML max r 2 )/ (1 + 
r 2 ), (—1 + / oc r 2 )/(l + r 2 ) and the values: 

(a) tx/M Lmin = 1.4096484 = D(l) 
M Lmin /M Lmax = 0.93 

= -0.93 

(b) t^/M Lmin = 1.5434 = D max 
M Lmin /M Lmax = 0.93 

foe = -0.93 

(c) tx/M Lmin = 1.4096484 = D(l) 

rrvin 

/M Lmax = 1/3 

ico = o 

Figs. 3 to 6 show /, a L , M E , M L , p E T, and p LS as functions of r, for these models. 



FIGS 3, 4, 5a, 5b, 5c, 6a, 6b, AND 6c GO HERE 



As expected, the areal radius has a minimum at r = and does not go singular, and both 
the Lorentzian and Euclidean 'densities' are well behaved at the transition. 

Model (a) has £ entirely inside the external event horizon R = 2Mi maX) as shown in Fig. 
5a, and it was found that there must be very little variation in M L (r) and f(r) (Fig. 3) to achieve 
this. The densities on each side of E are also only mildly varying, as shown in Fig. 6a. It could 
be thought of as a perturbation of a Kantowski-Sachs model, allowing the particle world lines to 
emerge beyond R = 2Ml briefly, before recollapsing back inside and encountering the signature 
change. Since this is true even for r — > oo, the particles do not fill the spacetime, and the model 
may be completed by matching to a vacuum exterior. This makes a very satisfactory model of 
signature change in the black hole topology. 

Model (b) differs only in having the largest possible value of t^/M Lmin . This results in 
ol, Me, Ml, Rs, pet, and p^s all having non-zero gradient at r = 0, meaning these quantities 
are discontinuous through the origin — see Figs. 5b and 6b. A feature of this model is that the 
bang time a^(r) has a maximum away from r = (Fig. 4), indicating the particle world lines self 
intersect somewhere in the time evolution of the Lorentzian part of the model [43, 42]. This 'shell 
crossing' is now a serious deficiency of the model, involving densities that diverge and go negative. 
Otherwise it is very similar to (a). 

Model (c) is complete since r — > oo covers the asymptotic regions of the model, and / — > 
means it is asymptotically flat. However the signature change surface passes out of R = 2ML m ax 
and extends to R = oo (Fig. 5c), which is not ideal, despite the nice density profile in Fig. 6c. 

Summary 

Within the Lemaitre-Tolman metric form, the wormhole topology is possible in both 
Lorentzian and Euclidean regions with and without matter (dust) present. The Euclidean re- 
gion bounces provided the mass function is positive (as defined in (44)), so negative mass models 
were not considered. Only constant t transition surfaces, which are orthogonal to the fluid flow, 
were considered. 



Vacuum to Vacuum signature transitions are equivalent to the uninteresting constant T 
Schwarzschild case. 

Signature transitions are possible inside the horizon if the density is non zero on both sides 
of the transition. Time symmetric Euclidean regions permit a second transition following the 
bounce, emerging into an expanding spacetime behind the past singularity. This is true both for 
the constant R transitions in uniform density Kantowski-Sachs type models, and the more general 
inhomogeneous case. The general case is particularly satisfactory as it models a finite cloud of 
dust. 

5. Conclusion 

We have succeeded in demonstrating the possibility that a change in the signature of space- 
time may occur in the late stages of black hole collapse, resulting in a Euclidean region which 
bounces and re-expands, passing through a second signature change to a new expanding Lorentzian 
space-time. The classical singularity at R = is thus avoided. Such transition surfaces necessarily 
have non- zero extrinsic curvature. 

The model of signature change employed here is strictly classical. Quantum cosmological 
questions, for example the relative probability of different sorts of transitions, have not been 
considered. We have based our notion of manifold continuity on the fulfillment of the Darmois 
type matching conditions, since they are invariant to the coordinates used, and no modifications 
are necessary to adapt them to surfaces of signature change. As discussed in [19], surface effects 
appear in the conservation laws, even when stronger conditions than Darmois' are imposed. 

Based on this approach, we have shown that signature transitions are possible in a spher- 
ically symmetric Lorentzian space-time, in both the Schwarzschild and Lemaitre-Tolman metric 
representations, though the ensuing Euclidean region might not be empty. Once the Israel iden- 
tities are adapted to signature change, continuous 'density' is no longer required. 

Within the Schwarzschild metric form, such a transition was possible on a constant T slice, 
but this can only span the outer region R > 2M. Conversely the constant R surface can be 
entirely inside the horizon, but does not lead to a bouncing Euclidean region. Thus these models 
are not satisfactory. 

A study of the geodesies in each region showed that the two Euclidean regions, R > 2M and 
R < 2M were in fact complete manifolds. It was found necessary to match geodesic 4-momenta, 
P M , at the signature change, in order that all geodesies could be continued. This naturally means 
P M and m 2 = \P^P^\ are discontinuous. This is consistent with the fact that the density can jump 
at a signature change. 

These results were generalised using constant t transitions in the Lemaitre-Tolman metric 
form. With suitable choices of the function /(r), this metric can reproduce the Kruskal-Szekeres 
topology of two sheets joined by a wormhole, but with non zero density. It also has a Kantowski- 
Sachs limit. It was found possible to have a signature change surface completely hidden inside 
the horizon R = 2M in the Lorenzian region, provided there was non-zero density in both the 
Lorentzian and Euclidean regions. In the Lorentzian region, the matter is of finite extent, and may 
be surrounded by vacuum. It was also possible for the Euclidean region to be time-symmetric, 
so that after the bounce, the matter expands through a second signature change into another 
Lorentzian region — a new universe. This makes a very satisfactory model of collapse, bounce and 
re-expansion of a mass concentration. Within the Lemaitre-Tolman form, a constant t signature 
change surface cannot be arbitrarily close to the Lorentzian singularity R = 0. One might expect 
such transitions to occur only a Planck time before R = 0, which would require us to consider a 
different equation of state in the Euclidean region. This may well relax the limits on rjj; and / 
that were found. 

It was found possible to hide the entire signature change surface inside the Lorentzian 



horizon R = 2M^, if the model is non-vacuum in the central regions, with a vacuum exterior. The 
matter is collapsing from not far outside the horizon, as may be expected for a collapsing compact 
object. The limit on / for a completely hidden surface implies (1) that all the infalling matter 
must be in a finite cloud, moving on tightly bound paths (R ma x < 4M^), surrounded by vacuum, 
and (2) that the black hole topology is required. This provides a classical bounce model of the 
kind we sought. A quantum cosmological analogue could be of interest in the context of Smolin's 
'natural history' of universes proposal. A treatment similar to that of Kerner and Martin [47, 48] 
could permit the non-zero extrinsic curvature that is required. 
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Tables 



Continuation 
Condition 


Signature Change Surface 


Const. R, R < 2M 


Const. T, R > 2M 


(i) u% = < 


Qe = qL, h 2 E = -h 2 L = 


q E = q L , h% + h 2 L = 2(1 -2M/R) 


(ii) u E = u L L 


h E = h L 


h E = h L 


(iii) P| = If 


2h E hl = (hl-h E )(2M/R-l), 
Qe = Ql, m E h E = m L h L 


2h 2 E h 2 L = (h 2 L + h%)(l-2M/R), 
Qe = Ql, m E h E = m L h L 



Table 1: Conditions resulting from the three continuation conditions at two types of signature 
change that retain a Schwarzschild metric form. 



Transition surface 


Const. R, R < 2M 


Const. T, R > 2M 


Allowed hi 


h 2 L >0>(l- 2M/R) 


h 2 L >{l- 2M/R) 


Allowed h E 


h'l < (2M/R - 1) 


h 2 E <(l- 2M/R) 



Table 2: Summary of allowed ranges of geodesic energy parameters on either side of the two 
signature change surfaces. 

Figure Captions 

Fig. 1. Diagram illustrating properties of combined Euclidean and Lorentzian geodesic 
paths. The curves give values of various parameters, as functions of hL for geodesies arriving at a 
constant T signature change surface at the point i?s = 4M(those curves starting at hi = l/y/2), 
and at constant R = 4M/5 surface (those curves starting at hi = 0). Vertical slices through the 
graph give the values of h L) h E) rriL/m E , Rj£ ax (where it exists) and R E rn or R E ax for individual 
geodesies. 

Fig. 2a. The function F(r) vs. r. 

Fig. 2b. The function D{rj) vs. rj. 

Fig. 3. f(r) vs. r for the three Lemaitre-Tolman signature change models, (a), (b) and 

(c). 

Fig. 4. a^(r) vs r, for models (a), (b) and (c). Note that (b) has its maximum away 
from r = 0, whereas (a) and (c) has it at r = 0. 

Fig. 5a. The run of R^(r), Ml(t) and M E {r) for model (a). 

Fig. 5b. Rn(r), Mt(r) and M E (r) for model (b). Note that none of these quantities are 
smooth through the origin r = 0. 

Fig. 5c. Rs(r) : M L (r) and M E {r) for model (c). 

Fig. 6a. The run of pi and p E for model (a). 

Fig. 6b. p L and p E for model (b). Note that neither quantity is smooth through the 
origin r = 0. 

Fig. 6c. p L and p E for model (c). 
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